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Abstract 

We reconsider the theory of magnetoresistance in hopping semiconductors. First, we have shown 
that the random potential of the background impurities affects significantly preexponential factor 
of the tunneling amplitude which becomes to be a short-range one in contrast to the long-range 
one for purely Coulomb hopping centers. This factor to some extent suppresses the negative 
interference magnetoresistance and can lead to its decrease with temperature decrease which is in 
agreement with earlier experimental observations. We have also extended the theoretical models 
of positive spin magnetoresistance, in particular, related to a presence of doubly occupied states 
(corresponding to the upper Hubbard band) to the case of acceptor states in 2D structures. We 
have shown that this mechanism can dominate over classical wave-shrinkage magnetoresistance at 
low temperatures. Our results are in semi-quantitative agreement with experimental data. 
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I. INTRODUCTION 



The problem of magnetoresistance in the hoping transport was addressed decades ago. In 
particular, an interest to this topic was related to important additional information provided 
by corresponding experiments (including estimates of the localization length). The most 
general and natural mechanism of positive magnetoresistance of orbital nature is related to 
shrinkage of the localized wave function by magnetic field; it was extensively reviewed in 
Then another important mechanism of orbital magnetoresistance was considered by Nguen, 
Shklovskii and Spivak (for the review see |2|). It is related to a presence of under-barrier 
scattering of hopping electrons by intermediate hopping sites and to interference between dif- 
ferent hopping trajectories. Note that for the effective interference the difference of lengths 
of different trajectories should not exceed the localization length which restrict the location 
of the trajectories to so-called "cigar region". A significance of this mechanism was empha- 
sized by the factor of exponentially-broad scatter of hopping probabilities corresponding to 
different "hopping resistors". As a result of "logarithmic averaging" over different config- 
urations the most important role is played by those interference patterns where the total 
hopping probability almost vanishes as a result of the destructive interference. The magnetic 
field suppresses the interference and thus the average effect is negative magnetoresistance 
which appears to be linear at weak magnetic field (although becoming quadratic at H — > 0). 
An important features of the approach discussed in 2] were as follows. First, the authors 
exploited an assumption of a presence of many intermediate scatterers. Second, following 
the theory [3], the authors assumed the preexponential factor to be equal to jj,/r where [i 
is scattering amplitude, r is a distance between the hopping site and the scattering center. 
The picture of interference magnetoresistance considered in [2] was very rich including a 
change of the sign of magnetoresistance, effects of spin glass etc. 

Somewhat later the problem was also discussed in [4], \^ where it was noticed that in 
realistic situations the number of intermediate scatterers is small and most probably equal 
to one or (in average) even less. Another important ingredient of the paper oj was a usage 
of wave functions typical for Coulomb centers which have not contained preexponential 
depending on r. In contrast to the "scattering states" of [2j which contained preexponential 
factors decaying with r, this situation can be specified as "strong scattering case". Note 
that, although in {5] the authors considered 2D hopping, they addressed to the case of 
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delta-doped layer and thus the asymptotic of the wave functions was similar to 3D. The 
important result of theory suggested in \^ was the following: the patterns of the negative 
magnetoresistance were almost universal predicting the maximum value of ~ 0.6 of the 
total resistance, and even the combination of the negative magnetoresistanse and positive 
wave-shrinkage magnetoresistance gave the maximum value of negative peak (with respect 
to average resistance) of around 40 percents. 

Unfortunately, these predictions for "strong scattering case" were not in a good agreement 
with experiment. First, in most of experimental studies the effect of negative magnetore- 
sistance have not exceeded 10 percents and typically was around several percents. Then, it 
was shown jf| that in 3D semiconductors the negative magnetoresistance is suppressed with 
a decrease of temperature after the crossover from Mott-type hopping (at higher tempera- 
tures) to Efros-Shklovskii hopping over the states within the Coulomb gap. In the paper 
6J we explained such a behavior as a result of a decrease of concentration of the scattering 
centers within the Coulomb gap. However our calculations were based on the assumption 
that the preexponential factor of the wave functions asymptotic corresponded to scattering 
states of {2 1 ("weak scattering case") rather than to hydrogen-like asymptotics exploited in 
j^]. Later 7] we have also demonstrated that to fit the experimental data one should also 
take into account spin mechanisms of magnetoresistance. The first one, considered in 8], is 
based on the fact that the intermediate scatterer should be occupied to produce a negative 
scattering amplitude. Thus the interference depends on the mutual orientation of the spin 
of the hopping electron and of the spin of scattering center. Without external magnetic field 
only one half of the configurations gives an interference. In magnetic field all localized spins 
are aligned which increases the role of interference and, correspondingly, leads to an increase 
of resistance. 

Another spin mechanism of positive magnetoresistance was first considered in Ref. |9j and 
then was studied in detail in [l0[ • It is related to a presence of doubly-occupied hopping sites 
(corresponding to the upper Hubbard band). Due to spin correlations on these sites requiring 
s-pairing of the spins (recall that we consider here electron rather than hole hopping) some 
hopping transitions are suppressed in magnetic field (like ones from single-occupied site to 
single-occupied site). 

As it was mentioned above, the incorporation of all of the relevant factors allowed us 
to reach a quantitative agreement between the theoretical model and experimental data. 
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However basing the scattering state asymptotic we exploited an assumption of correlated 
impurity configurations which had no solid theoretical prove. 

Another important request to the theory of hopping magnetoresistance was related to 2D 
hopping. As we have mentioned above, the theoretical model of [5] exploited 3D localized 
wave functions which do not hold for typical experiments for doped quantum wells where the 
wave functions have 2D character. Then, we should mention a new important experimental 



results 



ll| j 12( obtained for selectively-doped quantum well structures where both centers of 



the wells and centers of the barriers were doped ensuring a formation of the upper Hubbard 
band. These structures demonstrated a suppression of negative magnetoresistance with 
a decrease of temperature for the samples with higher degree of doping. Although we 
attempted to explain this behavior in a similar way as for 3D structures in jg], it hardly 
works because of an important difference between 2D and 3D physics. 

In what follows we will give a consistent description of magnetoresistance in both 3D and 
2D structures including different orbital and spin mechanisms. An important conclusion 
of ours is that in most occasions one deals with a "weak scattering case" rather than with 
"strong scattering case". If we are restricted to the lower Hubbard band, the decisive factor 
is related to the presence of charged centers outside of the "cigar region" not involved into 
interference. The random potential imposed by these centers restricts the extension of the 
hydrogen-like asymptotics of the scattering centers up to the distance to the closest charged 
center while outside this region the preexponential of the asymptotics appears to be similar to 
the one for the potential well case ("weak scattering limit"). For the case of the states within 
the upper Hubbard band an additional factor is related to the non-Coulombic potential of 
the scattering center which is also of a short-range character. The resulting picture of 
hopping magnetoresistance appears to be different from the one suggested in |5| (based on 
the pure Coulomb wave functions) and from the one of 3] (exploiting the assumption of 
large number of intermediate scatterers). We also emphasize a role of spin mechanisms of 
positive magnetoresistance which can dominate over wave- shrinkage magnetoresistance at 
low temperature. In this concern a special analysis is given to spin mechanisms for acceptor 
centers which have an important differences with respect to the earlier discussed case of 
donor impurities. 
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II. NEGATIVE MAGNETORESISTANCE IN 3D CASE. 



Let us consider negative magnetoresistance in 3D case. As it was mentioned above, an 
important ingredient to be included with respect to the previous studies is a random potential 
imposed by the intermediate charged centers (including both donors and acceptors). 

We shall start from a solution of a Schrodinger equation 

h 2 

+ U (t)V + £7(r)tf = Em. (1) 

2m 

Here Uq(t) is the potential of impurity (Uq = —a/r in the case of hydrogen-like impurity 
level) and U(r) is random potential that comes from the charged centers mentioned above, 
m is the electron mass in the conduction band (or a hole mass in valence band) and E is 
exact electron energy (we consider \E\ 3> U(r)). 

Because of the fact that typical hopping lengths are much larger then characteristic 
localization length a, we can solve |I| at r > a. Moreover, the typical hopping length 
rh appears to be much larger than typical distance between charged centers which can be 
roughly estimated as n -1//3 where n is a dopant concentration. Indeed, for 3D variable range 
hopping of the Mott type 

(T \ 1/4 

r h = £a, £=(^j , (2) 

where 

To * — 3 , (3) 
ga 3 

where g ~ nj Sb is the density of states, 8b being Bohr energy. Thus one obtains 

1/3 /2W^\ 1/4 
ThU 1 ~ I I (4) 

that is even for the Mott law r^n 1 / 3 >> 1. The more so it holds for the Coulomb gap regime 
where strongly exceeds the corresponding values for the Mott regime. 

Thus we are interested in asymptotics of the wave functions at distances much larger than 
n 1//3 which for moderately compensated material can be considered as the correlation length 
of the random potential imposed by the charged centers. If so, we can make an important 
conclusion. Namely, the random potential U is formed by the long-range Coulomb centers 
and in this sense the potential produced by the "parent" (for the considered wave function) 
impurity at distances larger than n -1 / 3 makes no difference with respect to potential pro- 
duced by other charged centers. In other words, for r > n^ 1 ^ 3 one should not discriminate 



between U and Uq and should assume that the resulting potential U has a spatial average 
equal to zero. 

Having in mind that r >> a we approach the problem of the asymptotics of the wave 
function by means of the WKB method. We introduce the function with an account of 
the normalization factor for the function \1/ as: ^ — (na 3 )^ 1 ^ 2 exp(— (j)/K). The Shroedinger 
equation in 3D case leads to 

(f)' 2 - h{4>" + 2(f)' /r) = 2m(C/(r) + \E\). (5) 

We will expand this equation into series with respect to k — > 0. In zero order we have 

(j>' = pm(U+\E\) (6) 

This order gives us the exponent. To get the pre-exponent factor we should use the first 
order of perturbation theory. Here we have in mind that the function (f)' Q at large r is actually 
a constant - the more so that its linear expansion in U(r) is averaged out. 
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Accordingly the expression for up to the first order is 

= J ^2m(U + \E\)dr + r/r mm . (8) 

where r min ~ n" 1 / 3 . And finally the wavefunction \l/ is 

1 



* = exp I ~n J v 2m{u + \ E ^ dr ) (9) 



Having in mind the considerations given above, we can average 

(l/h) J yj2m(U+\E\)dr = kr 

where 



k = (l/h) (^J2m(U+\E\) 

Note that k differs from k = \j2m\E\jh only in second order of U/\E\ (oc U 2 /E 2 ), as the 
mean value (U) is zero. Also we neglect U in the pre-exponent factor. 

So at distances from the scattering center larger than the correlation length of the ran- 
dom potential (assumed to be equal to average distance between the charged centers) the 



wave function asymptotics has a preexponential factor oc r 1 which agrees with the scheme 
exploited in [2J], I™, l7f for 3D hopping. 

n n 

Now, following approaches |2j, |6[ let us estimate the hopping probability between the 
sites 1 and 2 in a presence of intermediate "scattering center" with an account that the 
energies of the centers obey a relation \E 3 \ 3> \Ei\, \E 2 \ as 

P0C|J! + J 2 | 2 , Jl=I 12 , J 2 = -^B. (10) 

Here J\ and J 2 are hoping amplitudes related to direct and scattered path correspondingly. 
Note that the destructive interference (leading to negative magnetoresistance) implies that 
E3 < 0) which means that in the equilibrium the scattering site is occupied. 
The energy overlapping integrals are given as 

hi = Ej-^ exp(-r# 1 'a) (11) 

where we have assumed that > r m i n = n^ 1 ^ 3 ; Sb being the Bohr energy. Without 
a magnetic field this amplitudes are real. Though in the magnetic field their phases are 
different and hoping probability is 

P oc \J X + J 2 e i¥ f. (12) 

Here phase difference p is equal to tp = 27r$/$o> where $ is the magnetic flux through the 
surface bounded by hoping paths. $0 is the elementary magnetic flux. Accordingly, the 
interference magnetoresistance for the situation ip < 1 can be given as 



d s r 3 , (13) 



r(H) 
r(0) 

Here J = Ji + J 2 , g is the density of states and r 3 is scatterer position. Angle brackets 

corresponds to the ensemble average. We consider magnetoresistance to be determined over 

hops with small J, so we neglect the term J\/J in f|T3|) and get 

r(H) 
r(0) 



la '7W*-(f dE ^f la 



1+ w 



dh, . (14) 



J 2 (r :i ,E :i 

To obey Ji ~ J 2 one, first, should have r 12 ~ r 13 + r 23 with an accuracy of the order of 
localization length a. Then, having in mind the preexponential factors one notes that for 
small r m i n the only possibility to obey the relation is to have one of the distances, ri 3 or r 23 
to be small. We will assume that it holds for r 23 which is estimated as 

^23~^min^r (15) 
£/ 3 



Let us chose the surface at which J = and transform an integration over r% in a way 
d 3 r3 d 2 TZdR± where 1Z is the coordinate on the surface in question while R± is a coordinate 
along the normal to the surface where we assume that R± = corresponds to J = 0. In the 
lowest order in R± we have J = (dJ/dR±)R±. As it is seen, the integration of the logarithm 
term over R± gives 

ipJx 



dJ / dR\ 

Finally, the integration of the factor 

Ji 



dJ I dR± 

over d 2 lZ approximately gives a volume accessible for the site 3. Note that we have ri3+r23 < 
ri2 + a, thus the projection of r 2 3 to the plane normal to ri 2 should be less than (^a) 1 / 2 . 
As a result, the integration over the spatial coordinate r 3 gives 

r 2 f 2 

min B 



El 



'-tup (16) 



In its turn, the area of the interference loop (entering the estimate of ip) is 

ri2[r min —a) (17) 



' E 3 / 

Note that these estimates actually hold for all accessible values of r23 up to r23 ~ T\ij1. 
The final result depends on the behavior of g(e). For g = const (Mott-type hopping) the 
integration over E% is naturally controlled by the lower limit which accordingly to Eq. [T5l 
corresponds to the larger possible value of r 23 ~ ri 2 /2. In this case the r.h.s. of EqJTBI is 
oc r^ 2 where r h ~ r 12 . 

In contrast, for the Coulomb gap hopping the integration over E 3 is controlled by the 
upper level, Ec, corresponding to the edge of the Coulomb gap. In this case r.h.s. of EqJTBI 
is oc r h since the value of r 2 3 does not depend on r h . 

Now let us consider a combination of NMR with a positive magnetoresistance related to 
wave function shrinkage which can be estimated as 

ln^=(fV (18) 
p(0) \BJ y 1 

where 

B 2 = 19 

viae 1 
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Here a is a numerical parameter resulting from he percolation theory; for Mott type hopping 
a ~ 400 1] while for the Coulomb gap hopping different sources give a ~ 300 and a ~ 700. 
In its turn, NMR can be rewritten as 

where 

k = gM^Brminrh^cLa 1 ^ 2 Mott law 

3 5/2 

k = ^^- a E 3 r 2a 1/2 ES law. (21) 

Here k is the dielectric constant, r& is the typical hopping length for the states corresponding 
to the edge of the Coulomb gap while Ec is a width of the Coulomb gap. One sees that as 
a result we have minimum of resistance, 

H^=*B, ln^ = 4 (22) 

It is seen that the value of H min decreases with a temperature decrease irrespectively to the 
type of the variable range hopping. At the same time for samples corresponding to Mott law 
the absolute value of resistance in minimum increases with a temperature decrease while for 
the case of the Coulomb gap hopping it decreases with temperature decrease. 



III. NEGATIVE MAGNETORESISTANCE IN 2D. 

Let us now approach the problem of negative magnetoresistance in the 2D structure 
where impurity wave functions are quantized in the orthogonal to impurity plane direction. 
First we will consider the case when we deal only with single occupied or empty impurity 
centers (as it was done above for 3D case). 

Let us start with approximation of 2D impurity wave function in the r 3> a 
limit. Analogously to previous case we neglect Uo(r) and introduce function as ^ = 
(7ra 2 ) -1 / 2 exp(— (j)/Ti). The corresponding WKB equation is 

(f)' 2 - h{4>" + 4>'/r) = 2m(C/(r) + \E\). (23) 

Following the same procedure as was applied for 3D case we obtain 

* = exp (-i / y/2m(U+\E\)dr) (^) ^ . (24) 

9 



Analogously to 3D case this wave function is nearly equal to the potential well wave function 



oc exp(— kr)/y/r where k = ^2m\E\ + Uj /% which differs from k$ = y2m\E\/h only in 
the second order of U/E. 

Now let us consider negative magnetoresistance related to the interference contribution 
to the hopping probability. Following the same lines as for 3D case we obtain the equation 
similar to [TH except that the integration is over d 2 r 3 and the density of states g also corre- 
sponds to 2D. An important difference is related to the fact that now the value of J vanishes 
at 

?"23 = r min (£ B /E 3 ) 2 (25) 

With a similar transformation of the variables the integration of the logarithmic term 
over the coordinates gives 



&l\ (26) 

E 3 J 



while the effective loop area is 



/ ,£.2 \i/2 
ri 2 f r min y— J a J (27) 

Since in 2D in the Coulomb gap regime g oc e one notes that irrespectively to the hopping 
law the integration over E 3 is controlled by the lower possible values of E 3 leading finally to 
the estimates of r 23 ~ r%2- Thus one obtains 

k-2 = gM^Br l l ii n rh2a 1 ^ 2 a 1 ^ 2 Mott law 

k 2 = ^r£ 2 B r min r]! 2 2a 1 ' 2 a 1 l 2 ES law (28) 
e 4 

Thus, at is seen, for the situation considered above in 2D the only combination of interference 
NMR and wave- shrinkage PMR can not lead to a suppression of NMR with a decrease of 
temperature (at least for low temperature limit of linear NMR) since for both laws the 
temperature derivative of k stays to be negative. 

An important feature of the 2D quantum well structure is an easy possibility to have an 
occupation of the upper Hubbard band. Namely, if we dope not only the well regions, but 
also the barrier regions, the carriers from the barriers are captured by the wells and can 
form doubly occupied states. It was this situation which was realized in our experiments 
described in [12|. Since in these experiments we dealt with GaAs/AlGaAs structures 



of quantum wells with p-doping by Be, here we will also imply acceptor centers. 
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In our experiments the central regions of both wells and barriers were nearly equally 
doped by acceptor impurity Be. Thus the holes from the barriers have a possibility to 
occupy the second position for the acceptor in the wells forming the upper Hubbard band. 
However for the hole there was another possibility - to stay around its native acceptor in 
the barrier forming single-occupied center which we will denote as A . The corresponding 
scenario was first discussed in [13] . As a result, at the Fermi level we have centers with 
different occupation numbers - at least, A + (doubly occupied). A (single occupied) A 
(holes bound to the barrier acceptor) and A~ (empty barrier acceptor with no hole around). 

The possibility for the hole to form A + or A center depends on relation between the 
binding energies of these centers, Ub and Uf,. In particular, if [/& > Iff,, then all the barrier 
acceptors form A~ centers while all the acceptors in the well form A + centers. However for 
our experiments of the distance between the barrier acceptor and the interface between the 
barrier and well was not large and we expect Ub > Ub. In this case the probability to form 
A + center depends on the distance between the barrier acceptor and the closest acceptor in 
the well. Indeed, the formation of A + center profit from the interaction between A + center 
and A center 13] . 

Here we assume that some holes from the barrier are still coupled to their parent acceptors 
(A centers) and some are localized on the acceptors in the well (A + centers). According 
to charge conservation the number of A~ centers (that are free A centers) is equal to the 
number of A + centers. 

N{A + ) = N(A~). (29) 

In addition, we believe that there exists a random potential that overlap the energies of 
different types of centers. If the variances of A and A~ energies are equal, (129]) leads to 
equal densities of states for A and A~ at the Fermi level. For our purpose we assume that 
this densities of states are at least comparable. 

As for the negative magnetoresistance for the upper Hubbard band, it can be considered 
in the same way as for the lower Hubbard band discussed above. Note that the scattering 
potential strongly decays with distance Uq oc r -4 and thus the corresponding asymptotics 
of the wave functions are similar to the one given by Eqf54]but one should take r min = a. 
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IV. SPIN MECHANISMS OF MAGNETORESISTANCE FOR ACCEPTOR 
STATES 



We shall start from the mechanism of spin magnetoresistance first suggested in [8] which 
seems to be especially important for acceptor dopants. It is related to the fact that inter- 
ference can occur only if spin states of the final states for both tunneling paths coincide. 
For 3-cite configuration we discuss it means that the initial and intermediate centers should 
have the same spin projections (we remind that for destructive interference in question the 
energy of intermediate center should be negative, i.e. at the equilibrium this center should 
be occupied). For the case of acceptor states corresponding to the lower Hubbard band the 
corresponding configuration is in our case A — A — A~ where the role of intermediate center 
is played by the site A . Since the hole has spin 3/2, we have 4 projections of the spin and 
thus the probability for two sites A to have the same spin projections is P(H = 0) = 1/4. 
However in strong magnetic field the site spins are aligned and spin does not affect the 
(destructive) interference, that is in this case P(H — > oo) = 1. Thus an increase of the 
magnetic field leads to an enhancement of destructive interference which means positive 

n 

magnetoresistance which was noted in [8]. 

In a presence of the states representing the upper Hubbard band (in our case of A + and 
A states) the situation is somewhat more complicated. In particular, it is related to the 
fact that the spin structure of doubly occupied A + center is more complex than for single 



occupied site. In particular, the total spin of A + center is 2 (see 14j) and the possible spin 
states of A + center are the following: 

|J = 2,J. = -2> = i*V%-5*V9/» 



I j = 2, j, = o) = i*« *% 2 + \*[p% 2 - ±«V$ " I*-Vv* 
|j = 2,j. = i> = ±*«*%-i«V& 

|J = 2,J. = 2> = i»™»$-5*&»» 

where ty(}' 2 > - are the wave functions characterized by given spin projections of the two holes. 
Basing on these considerations one can show that for the destructive interference involving 
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purely the states of the upper Hubbard band, that is for configurations A + — A + — A , 
P{H = 0) ~ 1/4 while P(H -»• oo) = 1 

If the states of both of the Hubbard bands coexist at the Fermi level, it can be estimated 
that the average statistical factor P(H = 0) is still of the order of 1/4, although its value at 
strong fields, P(H — > oo) appears to be somewhat smaller than unity. 

At weak magnetic fields one expects a degree of spin alignment to be oc (figH) 2 /T 2 and 
thus the statistical factor is equal 

P(H)^P(H = 0) + a(^-) 2 (30) 



V T 

Here the coefficient a according to more detailed statistical calculations which we are going 
to present elsewhere can be estimated to be of the order of 1/2. 

Since P describes probability of the destructive interference, one concludes that at weak 
fields the positive magnetoresistance resulting from statistical factor P is quadratic in terms 
of magnetic field. It can be estimated as follows: 

AR sat 



R(0) \ T J 



R{0) 



(31) 



where and AR sat is the saturation value for the interference contribution to resistance with 
no account of spin degrees of freedom which is achieved when the phase (p in EqJT2l exceeds 

1 /2 

27T. As it can be estimated, the ratio AR sat /R(0) is oc r h for Mott hopping and oc r h for 
the Coulomb gap hopping. 

In its turn, it coexists with linear negative magnetoresistance (of orbital nature) which 
at relatively weak fields can be estimated as 

P(H = 0)AR sat - B 



H S at 

As it is known, if at the Fermi level the states of the lower and the upper Hubbard 
bands coexist, there also exists a specific spin mechanism of positive magnetoresistance 



first considered in [9j (and later discussed in detail in lOj) for n-type 3D structures. In 
such structures one deals with D° (occupied donors), D~ (doubly occupied donors) and D + 
(empty donors). Without external magnetic field the following configurations of hops are 
possible: D~ — > D°, D° — > D°, D° — > D + , D~ — > D + . In the magnetic field the spins of 
D° centers are polarized and thus the hops D° — > D° are forbidden (since in the final state 
of the second site corresponding to D~ the spins should be in opposite directions). In the 
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same way the transitions D~ — > D + are also suppressed. Thus the resistance increases as a 
result of application of external magnetic field. 

In our case of p-type structures the situation appears, again, more complicated due to 



3, 



more complex structure of A + centers. However, in general, the considerations given in 
Q still hold. Basin, on the calculations similar to given in [10] one obtains for weak field 
limit figH < T the following estimate: 

where C ~ 1/3, 

F = j^\ 2 (33) 
K9i + 9u) 2 

while gi, g u are the densities of states of the lower and upper Hubbard bands. Note that for 
the low concentration of dopants g u is controlled by the concentration of A + centers while 
gi - by the concentration of A~ centers and thus g\ = g u . At stronger magnetic fields when 
figH > T, the corresponding contribution to magnetoresistance still increases with magnetic 
field increase until figH reaches the value £T and then saturates jsj, jioj ]. 

One notes that at low enough temperatures the positive magnetoresistance of the spin 
nature suggested in [9[ can exceed the wave shrinkage magnetoresistance. At the same time 
this contribution at relatively weak fields when figH < T is expected to be comparable to 
the spin magnetoresistance resulting from interference term discussed above. Summarizing 
the both spin contributions to quadratic magnetoresistance we estimate the coefficient k 
resulting from the similar parametrization of the positive quadratic and linear negative 
magnetoresistance as was done above: 



k 2 = ^^wl^a 1 ^ ES law 



k 2 = gMEsr^inr^a 1 (3 Mott law 

,2 



e ^B> 

rp 3/2 1/2 a n 

P = P{H= 0)-^ r -^-^(CF + a^)-^ (34) 
gfiB en R(0) 

Thus, as it is seen, for the Mott case at T — >■ k oc T 1 / 3 while for the ES case it is oc T 1//4 . 

Note that in our calculations we assumed that the value of H min still corresponds to 
linear behavior of negative magnetoresistance which means that the magnetic flux through 
the interference area is much less than magnetic flux quantum <3> - The critical field H sat 
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corresponding to a crossover from the linear behavior to saturation of the negative magne- 
toresistance is given as 

H '" K i^fc (35) 

Correspondingly, if H min given by Eqf25] appears to be larger than H sat our calculations 
given above are invalid and one should compare positive magnetoresistance with saturated 
negative magnetoresistance rather than with linear negative magnetoresistance. One notes 
that in contrast to linear magnetoresistance which is proportional to the area of the in- 
terference loop for the saturation magnetoresistance this proportionality is omitted. As a 
result, as it was noted above, the temperature dependence of the saturation value of nega- 

1 /2 

tive magnetoresistance AR sat /R(0) results from factors oc r h for Mott hopping and oc r h 
for the Coulomb gap hopping. It is seen that the corresponding increase of the saturation 
magnetoresistance with temperature decrease is much weaker than increase of the positive 
magnetoresistance. Then, in the case H min > H sat it is the value of H sat which corresponds 
to minimal resistance since it separates a region of resistance decrease due to negative mag- 
netoresistance and resistance increase due to positive magnetoresistance. However at this 
situation spin magnetoresistance ( 13~T1) is also saturated so the temperature behavior of pos- 
itive magnetoresistance is related to (|18p and (or) to (|32p . 

In its turn let us consider the temperature behavior of the relation between R m in and 
H sat . According to Eqj22]and EqJM] 

j,2 r 5/2 

H min = 2g M E B r 1 J l la — -— oc T 7 ' & (Mott law) 

H ^ = ^ min a {g ^ )2CF r -£<xT (ES law) (36) 

At the same time H sat oc T 1 / 2 for Mott law and H s oc T 3//4 for ES law. Thus the ratio 
H min /H sat decreases with temperature decrease and this decrease is more pronounced for 
Mott law. 



V. DISCUSSION 



At FigJT] we present our experimental results from Ref. js] for 3D hopping concerning 
temperature behavior of magnetoresistance for regimes of Coulomb gap hopping ( Fig. [TJa) 
and of Mott-type hopping (Fig. CEJb). It is seen that these results are in in a qualitative 
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agreement with predictions of Sec. 2. In particular, the minimal value of resistance increases 
with temperature decrease for Mott type hopping and decreases for the Coulomb gap hop- 
ping. As it was noted in the Introduction, the agreement was strongly improved when we 
had taken into account more subtle spin effects [7|, however here we will not go into these 
details discussed earlier. 




FIG. 1: Temperature behavior of magnetoresistance for bulk CdTe crystals doped by donor 
impurities (CI). Fig. [TJa - the curves for the sample in the Coulomb gap regime, Fig(H b - for the 
sample in the Mott regime. 



At Figs. I2"f51 we present experimental results described in [12| for p-GaAs/AlGaAs:Be 
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multiple quantum well structures with different dopant concentration n. It is seen that for 
the sample with smaller concentration (Figj2]) the negative magnetoresistance is strongly 
enhanced with temperature decrease while for the sample with larger dopant concentration 
(FigfS]), in contrast, it is suppressed with a temperature decrease. 




FIG. 2: Temperature behavior of resistance for the structures of 10 GaAs wells of thickness 15 nm, 
separated by AlGaAs barriers with thickness 15 nm. The central parts of both wells and barriers 
were doped by p-type impurity Be with concentration 1 • 10 17 cm~3. 

To our opinion, the difference of magnetoresistance curves for samples with different n is 
related to the following fact. The sample with smaller concentration is far from the metal- 
insulator transition and the localization length is relatively small, a ~ lOnm. Thus, in a 
view of small n and small a the 3-cite approximation for interference contribution holds, 
rar/ l (ar m j n ) 1 / 2 < 1. For heavily doped sample a ~ 20nm and n ~ 10 12 cm -2 , correspondingly, 
nr/ l (ar m j n ) 1 / 2 > 1. As a result, the 3-cite approximation for this sample does not hold 
and the interference loop includes large number of scatterers. As it was noted above, the 
spin statistical factor for each additional site with non-zero spin at H = for acceptor 
impurities is ~ 1/4. Correspondingly, the interference contribution for loops involving many 
intermediate scatterers vanishes at H = 0. As a result, the linear contribution to negative 
magnetoresistance is, in any case, much smaller than for weakly doped samples. In contrast, 
the quadratic positive magneroresistance resulting from the statistical factor given by Eq. 
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FIG. 3: Temperature behavior of magnetoresistance for the structures similar to described at Fig. 
[2] but with concentration of Be 9 • 10 17 cm -3 . 

[3T1 strongly increases with temperature decrease, 

oc T- 7/3 (Mott law), oc T~ 9/4 (ES law) (37) 

In addition, we can expect that for the sample with large n H sat < H min , thus it is H sat 
which plays a role of H min . Due to weak temperature dependence of R sat the temperature 
behavior at the fields larger than H sat , that is corresponding to the minimum of p{H), is 
completely controlled by the spin PMR which gives p{H) oc T~ a with a > 2 . Indeed, an 
increase of resistance by a factor of 4-5 is observed at the fields larger than 0.3 T for the 
temperature variation from 0.9 to 0.4 K. 



VI. CONCLUSIONS 



To conclude, we have reconsidered existing theory of hopping magnetoresistance. We 
have shown that the random potential induced by the background impurities can affect the 
asymptotics of the localized states and, as a result, suppress to some extent the negative 
magnetoresistance related to interference effects. We have also generalized the theory for 
the case of acceptor states in 2D structures including the effects of the upper Hubbard band. 
The results obtained are in agreement with existing experimental data. In particular, we 
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explain the suppression of negative magnetoresistance with temperature decrease observed 
earlier for both 3D and 2D structures. 
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